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Abstract. In this paper the following three control systems for first-order 
time-delay plants are studied and compared: the feedback proportional-integral 
controller (PI), the Smith Predictor (SP) and a proposed variable structure 
consisting of two blocks. This structure acts as an open-loop proportional 
controller, after a setpoint change, and as a closed-loop integrating controller, 
when the error enters in a preset band. A chart, provided with the bor- 
derlines of the stability zone and with the curves of two design parameters, 
is implemented for each controller. The first parameter is the overshoot of 
the controlled variable, evaluated during a step change of the setpoint and 
made equal to a preset value. The second parameter, only for the PI and SP 
controllers, is the integral of the squared error (ISE), which must have the 
minimum allowable value. The ISE is also assumed as comparison index and 
the proposed controller appears as the best. 



1. Introduction 

First-order time-delay processes are exarained in this paper, since they are a 
simple but adequate approximation of most industrial plants. Moreover, controllers 
provided with two parameters are studied, since only systems, which can operate 
in a satisfying way when two performance indices are optimized, are taken into 
consideration. 

The structure, represented in Fig. 1 and provided with a proportional-integral 
(PI) controller, is the most used in time-delay systems thanks to its simplicity 
and to its robust performance in a wide range of operating conditions. Several 
tuning methods have been developed and their reliability and accuracy have been 
continuously increased. The first tuning rules, still prevalently implemented in the 
industrial systems since very easy to use, have been proposed by Zicgler-Nichols 
[1]. A better dynamical behaviour is assured by other methods, among them the 
performance criterion of Zhuan-Atherton [5] and the deadbeat response rule [5] 
deserve to be mentioned. Moreover, a new mathematical tool, consisting of the 
explicit solution of differential difference equations obtained by the method of steps 
[4], is recently available and therefore allows a more precise representation of the 
transient behaviour and an improved tuning of the control systems. 

In the past some attempts of finding new control implementations, suitable to 
eliminate some shortcomings occuring with the existing systems, have been made. 
The most significative is the Smith Predictor (SP) [5] (see Fig. 2), which, when 



Date: May 23, 2007. 

2000 Mathematics Subject Classification. 93C23; 34K35. 

Key words and phrases. Time-delay systems, integral of the squared error (ISE), overshoot, 
stability, variable structure controllers. 

1 



2 



GIANPASQUALE MARTELLI 



perfectly matched to the process, has the very useful property to act as a PI con- 
troller, subject to a delayed setpoint but operating in a feedback system without 
delay. 

In this paper a variable structure, depicted in Fig. 3 and based on already inves- 
tigated and well-known theories and procedures, is presented. It indeed makes use 
of both open-loop proportional and closed-loop integrating control, of the tracking 
function and of an adaptive technique updating the controller gain. 

The selected tuning indices are the overshoot of the controlled variable, named 
POy, and the integral of the squared error (ISE), both evaluated during a response 
to a step setpoint from the rated value to zero. It is also required an overshoot of 
the controller output, named PO^, lower than a preset value, in order to avoid the 
controller windup. The ISE, calculated for a definite time tg from the instant of 
the setpoint change, is then selected as comparison index. In order to reduce the 
computing complexity to an acceptable level, tg is assimied equal to seven times 
the process time delay, enough long to yield sound results. 

In Section 3 the new structure, described in Section 2, and the SP and PI con- 
trollers will be tuned and compared. All the information needed for their implemen- 
tation will be included in charts, having as coordinates the controller parameters 
and provided with the stability zone borderlines and with the curves of the perfor- 
mance indices, whose analytical expressions are detailed in Appendices A, B and 
C. Finally, Section 4 contains some conclusive remarks. 
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Figure 1 . PI controller 
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Figure 2. SP controller 



2. Proposed controller 

The proposed controller, depicted in Fig. 3, consists of two blocks, Ci and C2, 
selected respectively for the first and the second operation mode. It is necessary 
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to measure firstly the process gain, for example with an open- loop test, and to 
memorize it in the controller. The setpoint is continuously sampled and each value 
is compared with the previous ones, in order to check if its variation is inside a 
preset ztBg band (typically Bg = 2%). When the first operation mode is working, 
the C'2 output tracks the C'l output. 

When a setpoint change bigger than Bg occurs, the first operation mode, if 
not already working, is switched on, the controller acts as a open-loop proportional 
controller and the Ci output becomes equal to the value corresponding to the actual 
setpoint. When the error enters into the ±Bs band, the second operation mode is 
switched on and the controller acts as a closed- loop integrator. 

While the second operation mode is on and until a new setpoint variation re- 
quiring the first operation mode occurs, an approach of the controlled variable to 
the setpoint without oscillations outside the ±Bs band is necessary. This can hap- 
pen only if suitable values of the C2 parameter is chosen. When the setpoint is 
constant or stays inside a narrow preset band for a preset time, the ratio between 
the controlled variable and the controller output is memorized in the block Ci as 
updated process gain. 




Ci{s) 




Y(s) 



Figure 3. Proposed controller 



3. Main results 

A very useful implementation tool is represented by the tuning chart, provided 
with the stability zone borderline and with the curves of the performance indices. 
The tuning point is simply the intersection of two curves, each referred to a preset 
value of a performance index; the stability of the system under examination is 
assured, if this point lies in the stability region. These charts are plotted in Figs. 
4 (valid for tp — 0.55 as example), 5 and 6 and their coordinates are respectively 
{h]hi],[h]hitp] and [tp]hi]. 

During the response to the setpoint change from the rated value to zero, the 
variables r{t), u{t) and y{t), corresponding to i?(s), U{s) and Y{s), will be referred 
to their rated values. The SP tuning indices POy and POy are directly evaluated 
by finding the minimum of analytical functions in the range from t = to t = 00 
and are given respectively by (jB.7p and (jB.Sp . The index POy of the other two 
controllers is the opposite of the minimum negative of 701 values of y{tn) computed 
at equally spaced points; analogously POy is the minimum of v{tn). The used 
expressions are (|A.3P and (|A.4|) for the PI and (|C.3|1 or (|C.4|1 and (|C.5|1 for the 
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proposed controller (1 <= n <= 7); the corresponding curves Ty {POy = 0.0105) 
and Vv {POv — 0.1) are plotted in the tuning charts (see Figs. 4 and 6). The ISE 
of each controller is evaluated with the triangle rule for the same time (from t = 
to t = ts = 7), according to CO|) . (|R9)) and ((CJ)) . 

3.1. Transfer functions. The transfer functions of the blocks included in Figs. 1, 
2 and 3 are given by 

(3.1) p(,)^p„(,)e-i = ^,^e-^^ 

1 + sTp 

(3.2) C(s) =Kp + ^ 

(3.3) C2(s) = — 

s 

where K represents the plant steady-state gain, Tp the plant time constant, L the 
plant time delay and Kp, Ki the controller parameters. It is convenient, in order 
to get equations independent of the real values of the parameters, to introduce the 
normalized time referred to the plant time delay L {<j — s L) and the dimensionless 
parameters tp = Tp/L, h — K Kp and hi = K KiL. 

Considering a matched SP (Mo(ct) — Po{cr) and M{a) — P{a)), the relationships 
between the controlled variable Y{a), the setpoint R{(j) and the controller output 
U{a), based on the transfer functions from p.ip to p.3p . are the following: 

(1) PI controller 

(3.4) Y{a) = Cia)P{ama) Yia)) = e-'^{R{a) Y{a)) 

(t(1 + atp) 

(3.5) U{a) = Y{a) = e^Y{a) 



(2) SP controller 

I + C(c7jPo[a) a-^tp + all + hj 



(3.7) U{a) ^ Y{a) = Y{a) 

(3) proposed controller - first operation mode 

(3.8) Y{a) = P{a)U{a) = U{a) 

(4) proposed controller - second operation mode 

(3.9) Y{a) = C2{cT)P{ama) - Yia)) = ^' e-'^(i?(<7) - y(a)) 

(t(1 + crip) 



(3.10) 
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3.2. Stability. The stabihty of the SP controllers is theoretically assured for 1 + 
h > 0, which always holds since h > 0. Thr instability may really happen if the 
process model transfer function is not correct, but in this paper this is excluded. 
The following, proved in [6^ and confirmed in [7j, hold for the PI {h 0) and the 
proposed {h — 0) controllers: 

(1) The parameter h must be included in the range from zero to given by 

hu = - COs(Za) + tpZaSin{Za) 

where Za is the first positive solution of 



tan{za) 



1 + ip 



(2) The parameters hi, for a given h, must satisfy the following inequalities 
5r{zi) < anddr{z2) > and hi > 

where 



Sr{z) — — [hi — z sin{z) — tpZ^cos{z) 



1 
L 

and Zi and Z2 are the first two positive roots of 



h + cos(z) — tpZ sin{z) — 

It is also useful to take into consideration the stability phase margin (PM), whose 
expression is 

„ N hi + z?ht„ 

tan(zb + PAf) = 

Zb[h - hitp) 

where Zb is the solution of 

h"^ 

h^ + -^ = l + tlzl. 

The curves corresponding to the stability region borderline and to PAI — 30°, 
PM — 45° and PM — 60°, named respectively Tg, Fpi, and F^s, are plotted 
in the tuning charts (see Figs. 4 and 6). 



3.3. Saturation. One of the nonlinear issues of any controller is the actuator sat- 
uration, since from the practical point of view the actuator transfer function is 
bonded and the controller output may be out of the rated band. This windup is 
avoided, in the simplest way, by switching off the controller integrating action when 
the actuator saturates, but this mode introduces some control discontinuities. A 
more sophisticated solution is given by the addition of the anti-windup scheme, 
depicted in Fig. 10.50 of [8]. An other simple solution without discontinuities, 
adopted in this paper, consists in setting the rated value of the controller output 
lower than the extreme point of the linear part of the actuator transfer function 
and in tuning the controller in such a way that this point is never reached. 
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3.4. Steadiness. When the proposed controUer is switched from the first to the 
second operation mode, the controUed variable, approaching the setpoint, should 
enter into and remain inside the preset ±_Bs band (typically Bg = 0.02) until a new 
setpoint variation requiring the first operation mode occurs. This happens if the 
parameter hi is enough low and precisely if the tuning point is below the curve Tf,, 
corresponding to POb — Bs, where POb is the maximum absolute of 70 f values of 
y{tn), computed at equally spaced points with (|C.3|1 or (|C.4p (see fig. 6). 

3.5. Tuning. 

(f) PI controller 

The tuning is performed according to the deadbeat response rule, presented 
in [3] and summarized as follows: the ISE must be minimum and the 
controlled variable overshoot POy must be equal to a preset value {POy = 
0.0105), provided that the controller output overshoot PO^ is not higher 
than a fixed value {POy < 0.10). The tuning point B, having the minimum 
ISE value and lying on Ty, is indicated in Fig. 4. 

(2) SP controller 

It is tuned with the same criterion as the PI controller. Since the POy and 
ISE values on each point, lying on Fy, increases and decreases respectively 
when the coordinate h increases, the intersection of Ty and F„, named B, 
must be assumed as tuning point. The damping borderline F^, correspond- 
ing to 6 = according to (jB.3p . is also plotted in Fig. 5; the underdamped 
setpoint response happens for points above F^. 

(3) proposed controller 

This controller has one parameter (hi) and, therefore, it is sufficient to 
require a controlled variable overshoot equal to a preset value {POy = 
0.0105). The point, lying on the related curve F^ in Fig. 6 in correspondence 
of the coordinate tp, can be considered as tuning point, since below both 
the controller output curve F„ {POy = 0.0105) and the steadiness curve F{, 
{POb = Bs = 0.02). 

The tuning parameters h and hi and the performance indices ISE, related to some 
values of tp from tp = 0.1 to tp = 10, are gathered in Table [1] (as per Tables 1. 
and 2. of [3] for the PI controller). It is plain that the proposed controller has the 
minimum values of the ISE for all tp and therefore it can be elected the best. 

4. Conclusions 

In this paper a new controller is presented. Its main characteristics are: 

(1) during the switching from the first to the second operation mode, it is 
bumpless since the controller outputs of both blocks have the same value 

(2) it is adaptive, since the process gain is continuously monitored and accord- 
ingly updated, when a discrepancy appears between model and process 
values 

(3) the allowable values of the controller parameters, which avoid the actuator 
saturation, are easily deduced from the tuning chart 

(4) the process mismatch can be easily taken into consideration in the tuning 
chart, since the process parameters are included in their coordinates and 
suitable margins can be previously selected 
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B : tuning point; Ts : stability borderline 

Tpi : PM = 30°;Tp2 : PM = 45°;Tp3 :PM = 60° 

Ty : POy = 0.0105; : POy = 0.1 



Figure 4. PI controller tuning chart for tp = 0.55 
Table 1 . Parameters and ISE of the controllers 
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Figure 5. Smith Predictor controller tuning chart 



(5) it has zero steady-state error during the response to a step setpoint also for 
an integrating process 

(6) the lowest values of the ISE can be achieved 

(7) it can be easily implemented in the modern digital controllers. 

The proposed controller may replace the PI controller, especially when the system 
is manually controlled, significative setpoint changes are needed and lowest ISE 
values are required. 



Appendix A. Response to a step setpoint of the PI controller 

During a response to a step setpoint from r = 1 to r = applied to the closed- 
loop system in steady condition at i = 0, the differential equation, related to the 
controlled variable y and deduced from (|3.4p . is given by 

dy{t) (fy{t) 
(^p^ -y-^ dt ^ dt"^ 

^ih, + h^){r{t^l)-y{t-l)) fort>=0 
dt 

subject to an initial condition of the form 

(A.2) 2/W = 1 for -l<=t<=0 
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Figure 6. Proposed controller tuning chart 



The analytical solution of (|A.1[) . obtained with the method of steps in [1], consists 
of a set of functions t/„ (t„ ) , each valid for n ~ 1 < t < n and given by 

yiih) - 1 

/ \ o\ i=n~l j=n-2 

1=0 j=0 

where t„ has the temporal origin at i = n — 1. In Example no. 2 of |4j there are 
also the recursive expressions suitable for the evaluation of the coefficients An^i and 

Considering t„ = ^n-fi, from (j3.5p one obtains for Vn {tn) = K Un{tn) 

(A. 4) Vn(tn) = Vn+l{tn+l) + — 

atn+1 

Finally, the integral of the squared error is given by 

ISE = + ]^ 0.01 {vM? - \ 0-01 {ytX^)f 

(A. 5) n=t^ r=100 

+ 0-01E E (y"(o-oi^))' 

n— 1 r— 1 

where is — 7. 
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Appendix B. Response to a step setpoint of the SP controller 

During a response to a step setpoint from r = 1 to r = applied to the closed- 
loop system in steady condition ai t = 0, the equation, related to the controlled 
variable y and deduced from l|3.6p . and its underdamped solution are given by 



(B.l) 

(B.2) 
where 



tp^-^ +{l + h)^ +h,y{t)^h, forO<t<l 



t,^^ +{l + h)^ +h,yit)^0 fort>l 



y{t) = 1 for < < < 1 

y{t) = e-''^*-^\cos{h(t-l))+'^ sin{h{t-\))) for t>l 
^_ 0.5(1 + 

tp 

(B.3) h^^-^ {Ah,t^~{l + hYf-' 

bp 

From ((3Jl) and (|R2)) one obtains for v{t) = Ku{t) 
v(t)^y{t)+tp'^ 

= {cos{b {t-l))- -^ + ^P^^ + ^^) _ 1))) 

The derivatives of y{t) and v[t) with respect to t are given respectively by 
(B.5) =_e-(*-i)^,,,(6(t_l)) 

(B.6) ^ =e-'^'^'-^\a' + h''){-Ucos{h{t-l))+'^^^^ sin{h{t^l))) 
at h 

Introducing the solution h{t — 1) of (|B.5[) . equated to zero, in (|B.2p . we obtain for 
the overshoot POy 

(B.7) POy = -e-'"'/'' 

Analogously introducing the solution h{t — 1) of (|B.6p . equated to zero, in (jB.4p . 
we obtain for the overshoot POv 

(B.8) PO, = ~e-^'''\\ - 2 flip + t2(a2 + feS^^o.s 

where 4> is the first positive root of tan{(f>) — btp/{atp — 1). Finally, the integral of 
the squared error is given by 

(B.9) ISE^l + ISEa{ts)-ISEa{l) 

where ts = 7 and 

+ a^(a^ — 3 h^)cos {2b{t- I)) + ab{-3 a"^ + b^)sin{2b{t - 1))) 
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Appendix C. Response to a step setpoint of the proposed controller 

During a response to a step setpoint from r = 1 to r = applied to the closed- 
loop system in steady condition at i = 0, the following holds: 

(1) First operation mode 
From (|3.8p one obtains 

ya{t) = l forO<=t<=l 

^^'^^ ybW - e(-*+i)/*^ fort>=l 

since yb{t) is the solution of 

2/h(i)+ip-^ =0 

subject to the initial condition yb(l) = 2/a(l) = 1- Considering the identity 
Vbi^ + tq) = Bg, the time tq required by j/f, to reach Bs is given by tq = 
tpln{l/B,). 

(2) Second operation mode 
From (13.911 one obtains 



(C.2) 



dy{t) d^y{t) 



= {K + h- ){r{t - \) - y{t - I)) fort>=l + tq 



dt " " dt^ 

d_ 

dt 

subject to the initial condition 

(a) tg <1 

y{t) = ya{t) = 1 for tq<t<l 
y{t) = yb{t) = e(-*+i)/*- for K t < 1 + tq 

(b) tq > 1 

yit) = yb{t) = e'-'+i'/*" for tq < t < 1 + tq 

The analytical solution of (jC.2p . obtained with the method of steps in 
[3], consists of a set of functions yn{tn), each vahd for n >= 1 and given by 

(a) tq<l 

i—n j—n—1 
ynAtn) =^An,l,^ti^+e'^"^*'' -^"^J*" < t^ < I - tq 

i—n—1 j=n 
ynAtn) = ^".2,itji + e-*"/*" Bn^.jti f Or 1 - < t„ < 1 

(b) tq > 1 

i—n—1 j—n 
(C.4) yn{tn)= Y "^nA+e-'-^'-Y^^-l^n f OT < i„ < 1 

where f„ has the temporal origin aX t = tq + n. The initial conditions are 
given by 

(a) tq < 1: yoAto) = 1 and yoMto) = e-*"/*^ et-'^+i)/*" 

(b) > 1: yo{to) = e"*"/*" e(-*'+i)/*p 



(C.3) 
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In Example no.l of [4j there are also the recursive expressions suitable for 
the evaluation of the coefficients An,i, Bn,j, An,is, Bns.j, An,2,i and Bn,2,j- 
Considering i„ = from (j3.I0p one obtains for = K Un{tn) 

(C.5) Vnitn) = +tp — 

atn+i 

Finally, the integral of the squared error is given by 

ISE ==1 + ^ (1 - + /S'£:6 /orO <tq<ts-l 

(C.6) 2 

ISE =1 + I (1 - e-2(*-i)/*P) forts -Ktq 

where tg = 7 and ISEi,, evaluated from t = 1 + tq to t = ts with an expression 
similar to (jA.Sp . is the contribute of the second operation mode. 
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